Acoustic attenuation in glasses and its relation with the boson peak 
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A theory for the vibrational dynamics in disordered solids [W. Schirmacher, Europhys. Lett. 
73, 892 (2006)], based on the random spatial variation of the shear modulus, has been applied to 
determine the wavevector (k) dependence of the Brillouin peak position (ilk) and width (Tk), as 
well as the density of vibrational states (g(u)), in disordered systems. As a result, we give a firm 
theoretical ground to the ubiquitous k 2 dependence of Tfe observed in glasses. Moreover, we derive 
a quantitative relation between the excess of the density of states (the boson peak) and I\, two 
quantities that were not considered related before. The successful comparison of this relation with 
the outcome of experiments and numerical simulations gives further support to the theory. 

PACS numbers: 65.60.+a 
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The most striking differences between glasses and crys- 
tals, at a macroscopic level, concern the thermal proper- 
ties. At few tens of Kelvin, the specific heat of glasses 
exhibits an excess over the Debye expectation, which, 
in a C(T)/T 3 vs. T plot, appears as a characteristic 
maximum. Similarly, at low-T the thermal conductivity 
increases as T 2 with increasing T and, in the same tem- 
perature range where C(T)/T 3 has a maximum, exhibits 
a plateau. While the first observation can be ascribed to 
the presence of an excess of states over the Debye den- 
sity of states 3d(w) (the "boson peak", i. e. the peak 



observed in g(u>)/u> 



w), the second one was until 



recently not understood at all. Both anomalies appear 
to be strongly affected by the characteristics of the nor- 
mal modes of vibrations in the THz frequency region. 
The nature of these modes has been the subject-matter 
of a very intense and controversial debate in the liter- 
ature since many decades. The issue has been recently 
revitalized thanks to new neutro_n, _X^ay_j _and other in- 
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The boson peak shows up in a frequency range where 
the broadening of the acoustic excitations becomes of the 
order of magnitude of the resonance frequency, thus in- 
dicating a possible mode localization. This observation 
lead different authors to hypothesize a relation between 
the position of the boson peak and the existence of local- 
ized vibrations. Acoustic waves that become Anderson- 
localized, indeed, could produce the plateau in the ther- 
mal conductivity. Following this idea, investigations of 
(Anderson-)localization properties of waves in disordered 
systems based on simulations model calculations 

[l^ . [2fjj and field-theoretical techniques [28| have shown 
that Anderson-localized states in disordered media do ac- 
tually occur, but in a much higher frequency range (near 
the upper band edge) than the boson peak frequency. 

So the question is still open: what is the very nature of 
the states near and above the boson peak frequency? As 
these states are neither really propagating nor localized, 
Feldman and coworkers [12j suggested to call them "dif- 
fusons" (they behave like diffusing light in milky glass). 



In this regime, however, the resonance frequency ilk still 
exhibits a linear dispersion with the wavevector fc, al- 
beit with an apparent sound velocity being somewhat 
larger than the ultrasonic speed More interesting, 

in this frequency range the width Tk of the excitations, 
a quantity proportional to the sound attenuation coeffi- 
cient, shows a k 2 dependence in most of the investigated 
materials [6] , whereas at lower frequencies the width ap- 
pears to have a stronger k dependence @ . These findings 
are still awaiting a theoretical explanation. 

In this letter we exploit a recently developed theory of 
vibrational excitations in disordered elastic media [2£| 
and apply it to the calculation of the dynamic struc- 
ture factor. This theory is based on the model assump- 
tion that the disorder leads to microscopic random spa- 
tial fluctuations of the transverse elastic constant (shear 
modul us) [301) . As in similar, more schematic approaches 
[2Cll23l ,l2a] the excess DOS has been shown to arise from a 
band of disorder-induced irregular vibrational states, the 
onset of which approaches lower frequencies as the disor- 
der is increased. Within this framework it was possible 
to formulate a theory for the energy diffusivity, which 
gave the first explanation of the plateau of the tempera- 
ture dependent thermal conductivity and its relation to 
the excess DOS. When applied to the calculation of the 
shape of the Brillouin resonance, this theory explains i) 
the ubiquitous k 2 dependence of Tk and ii) the observed 
increase of the apparent sound velocity with frequency 
[3l| . Moreover the theory predicts Hi) that the excess of 
vibrational states, if properly normalized with the Debye 
DOS, takes almost a universal value and iv) the existence 
of a quantitative relationship between the excess over the 
Debye DOS and the width of the Brillouin line. The lat- 
ter, previously unexpected, relation is shown to agree 
with the findings of experimental inelastic scattering and 
computer simulation experiments. 

We now shortly summarize the theory. We 
consider an elastic medium with a mass density 
mo, shear modulus G, bulk modulus K = A + |G 
(A = longitudinal Lame constant). These elastic con- 
stants are related to the longitudinal and transverse 
local sound velocities as c 2 -^=G/mo; c\ = (K + |G)/mo 
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FIG. 1: a) Reduced DOS (g(oj)/g D (oj) for e = 1 -70/70 
= 0.0001, 0.001, 0.01 and 0.1 (from left to right) and for 
(cl/ct) 2 = 2 (full lines) and (cl/ct) 2 = 6 (dots), b) Reduced 
excess DOS Ag(uj)/gz>{uj) for the same parameters [33|. 



= (A + 2G)/mo. We now assume that the shear modu- 
lus (and not A, which is set to a constant value Ao) ex- 
hibits a random spatial variation: G(t)—Gq[1 + AG(r)]. 
The random function AG(r) is supposed to be Gaus- 
sian distributed around its average Go with a variance 
oc 7g- The parameter 70 describes the "degree of disor- 
der" of the system. We emphasize that this phenomeno- 
logical model may be adequate both for a topologically 
disordered system (glass) or a disordered crystal. Us- 
ing standard field-theoretic techniques [HI, [H, IHJ 
the self- consistent Born approximation (SCBA) for the 
(complex) self energy function S(w)=E'(o;) + iE"(oj) has 
been derived (29|, resulting in the set of equations: 



E(w) 



: lG^k<k D [XL(k,uj) 

r 2r . .2 , ,2/2 



XL(k,uj) = k 2 [-cj 2 + k\c 2 Lfl - 2EH)]- 1 
X T(k,uj) = k 2 [-uj 2 + k 2 (c 2 - EM)]- 1 . 



(1) 



Here \l (xt) are the longitudinal and transverse dy- 
namic susceptibilities, resp., and cl.o, cu,o are the sound 
velocities of a system without disorder. The sum over 
k indicates integration up to the Debye cutoff, kn = 
[6ir 2 N/V] 3, according to the identification Efc<fc D — > 
(3/ftf,) f Q feD fc 2 d/c. The SCBA is the simplest form of an 
effective-medium theory for the disorder [lj|, [2(| [2^1 , in 
which the disorder effects lead to a frequency dependent 
modification of the sound velocities ("complex acoustic 
indices of refraction") by the function E(u>). The real 
part E'(u;) causes a renormalization and dispersion of 
the sound velocities according to cl(cj) 2 = c\ M — 2E'(w) 
and ct(uj) 2 = ct,o — E'(w), whereas its imaginary part 
describes the sound attenuation (see below). Apart 
from the natural length and frequency scales k~^ and 
= C£>kr> [33j there are only two non-trivial parame- 
ters in this theory, namely the degree of disorder and 
the ratio of the renormalized sound velocities Cl/ct [HI]. 



FIG. 2: Excess DOS g(u) — <7d(w) for the same e values as in 
Fig. [T] and for (cl/ct) 2 ~6 (full lines) together with the Bril- 
louin linewidth /(cl, ct)T(^)/^d= 2oj/(cl, ct)^" {uj)/lodc 2 l 
(dotted lines.) 

Insert: Onset frequency (full line), boson peak frequency 
(dashed line) and Ioffe-Regel frequency ujir — CL(io)k/-K 
(dash-dotted line) as a function of the parameter e= 1— 7g/7g- 



From Eq. (fT|), the DOS can be calculated as: 

g(oj) = (2w/37r)E fe<feD pIm{x i (fe,a;) + 2xr(A; ! a;)} (2) 

As in other theories of a harmonic solid with quenched 
disorder [2(| [H, [26|] the system becomes unstable if the 
disorder exceeds a critical value 7 G , which — in the 
present case — slightly depends on the ratio of sound 
velocities and ranges, for example, from 7^=0.1666 for 
c 2 /c 2 i=2 to 7^=0.227 for c\jc\=%. It has been demon- 



strated [2fJ, |23J, |25|, I2J, |29(, that the quenched disorder 
produces an excess over the Debye DOS, the offset of 
which approaches lower and lower frequencies as 7g^7g- 
In Fig. [T]we report the reduced DOS g(w) / gr)(oj) (i. e. 
the usual boson peak representation) and the reduced ex- 
cess DOS [g(u>)/ gr)(u>)]— 1 against uj/ujd for several values 
of the separation parameter e=l— 7g/7 g . It is remark- 
able that the value of the normalized excess does not ex- 
ceed the value 1, which it takes for small enough values 
of e. It is also remarkable that the excess essentially does 
not depend on the ratio of the sound velocities. More- 
over, the excess is significantly different from zero only 
above a certain frequency threshold (onset frequency), 
below which the DOS basically coincides with the Debye 
expectation. The excess DOS turns out to vanish as w 4 
for to — > 0. 

We now turn the attention to the dynamical structure 
factor, S(k,uj), which can be measured by inelastic neu- 
tron, X-ray and light scattering: 

S(k,tj) = i(n(w) + l)Im{Xi(fc,w)} = -("(">) + 1) 
x2fc 4 E"H {[k 2 c L {io) 2 - lo 2 ] 2 + 4fc 4 E"(cj) 2 } -1 
! r , s ^k 2 fc 2 E'»/co 

7T 2ld [cl{U))K — U)\ Z + [k Z L,"(UJ)/Lj\ Z 
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MD-LJ 


39.5 


1320 [16] 


540 [34] 


2.1 10"^ [34] 


0.120 [16] 


1.14 


610 


6.99 


8.8 10" 


-a 


2.4 


0.099 


0.103 
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MD-Si0 2 


14.1 


5940 [17] 


3570 [17] 


3.7 10" 5 [IT] 

5.3 10 -5 [35] 

1.4 10" 4 [3JJ 
9.9 10" 5 [39] 

3.5 10" 4 [41] 
2.2 10 -3 [43] 


0.045 [17] 


1.61 


3950 


63.7 


1.2 10" 


-5 


3.1 


0.00055 


0.00060 
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Si0 2 


15.1 


5960 [9] 


3750 [9J 


0.042 [4] 


1.58 


4130 


65.1 


1.1 10" 


-5 


4.9 


0.00101 


0.00053 
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Ge0 2 


16.1 


3650 [36J 


2150 [36] 


0.034 [37] 


1.54 


2380 


36.8 


6.0 10" 


-5 


2.3 


0.00237 


0.00187 
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Glycerol 


8.2 


3600 [38] 


1870 [38] 


0.031 [40] 


1.93 


2090 


40.4 


4.5 10" 


-5 


2.2 


0.00177 


0.00194 
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10.7 


2940 [41] 


1370 [41] 


0.036 [40] 


1.77 


1540 


27.3 


1.5 10" 


-4 


2.4 


0.00445 


0.00523 
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Selenium 


30.0 


1800 [42J 


895 [42] 


0.070 [44] 


1.25 


1000 


12.6 


1.5 10" 


-3 


1.4 


0.00935 


0.02170 


80 



TABLE I: Data taken from the literature together with quantities derived from them: atomic volume V/N, longitudinal and 
transverse sound velocities cl,t; coefficients A, Ad of the w 2 dependence of the measured and Debye DOS; coefficient B entering 
into the frequency dependence of the Brillouin width, Y{lo)—Bu 2 ; Debye wavenumber kD = [Qi^ 2 N /V] 1 ^ 3 , Debye sound velocity 
c D [H[ and frequency uj d = c D k D ; p (=4/4 D ), R exp =(A ~ A D )/B, R Th =f(c L ,c T )/u 2 D , Err=\R Ex p-RTh\/[(REx P + R T h)/2]. 



From this equation we can read off the width of the Bril- 
louin line (full width at half maximum) as: 

r(w) « 2fc 2 £"H/^ a — ^ wE"(w) (4) 

As the "disorder function" enters also into the 

DOS we find from © for < c l ,tM 2 the following 

approximate relation: 

uj D Ag(u) = w D [g(uj) - 9d{^)\ = f{cL,c T ) (5) 

LU D 

with f(c L ,c T ) = -( CD /c L ) 2 [l + (cl/ct) 4 } (6) 

To check the validity of the approximation made in de- 
riving Eq. ([5]), in Fig. [2] we compare the excess DOS 
with the Brillouin linewidth, multiplied with the factor 
/(cl,ct). We see that within our theory the approxi- 
mate relationship given in Eq. ([5]) holds, particularly in 
the regime where both quantities follow approximately 
an ui 2 dependence. 

Equation {5| represents the main result of the present 
work. It tells us that the Brillouin width T(uj) is propor- 
tional to lu 2 (and thus to k 2 recalling the slightly chang- 
ing value of Cl(u))) in the whole frequency region covered 
by the boson peak, and that it turns to a cj 4 behav- 
ior only below the onset frequency. The onset frequency, 
that must not be confused with the boson peak frequency 
ujbp (see insert of Fig. lies well below lobp, it is lo- 
cated at the low frequency edge of the boson peak (i. e. 
at that frequency where the excess DOS vanishes) and it 
tends to zero as e— >0. Equation ([5]) also gives a quan- 
titative prediction: it establishes a relation between the 
absolute value of the excess DOS and the Brillouin peak 
width. The material dependent physical quantities en- 
tering into this relation are only macroscopic properties, 
as sound velocities and density, that can be determined 
independently. 

In order to test the theoretical prediction manifested 
in Equation ([5]) we have collected a number of experi- 
mental and numerical simulation data of the vibrational 
DOS and the Brillouin line width, as well as sound veloc- 
ity data (Table I). A complete set of data, to our knowl- 
edge, only exists in the case of the numerical simulation 



of Lennard- Jones particles and for a model of Silica glass, 
and for the experimental cases of Silica, Germania, Sele- 
nium, Glycerol and o-Tcrphenyl. 

The first test is that of Fig. 1 a) from which we expect 
that the quantity p = [<7(w)/<7i)(w)] max = A/ Ad (where 
A, Ad are the coefficients of the measured and Debye u> 2 
law in the boson peak regime) is equal to 2. A systematic 
variation upwards is observed, but the overall set of data 
is not far from this value [45|. Exceptions are the case 
of silica, where the presence of specific localized modes 
(tetrahedra rotation) are well documented, and that of 
Selenium, where the ratio is low, most likely because of 
the boson peak is not well developed, indicating a large 
value of e. 

The second test involves the quantity 
R exp =Ag(uj)/r(Lu) = (A - A D )/B, which is the ra- 
tio of the excess DOS to the Brillouin width as 
determined by the experiments, while Rth is the same 
ratio as determined from the theoretical prediction, i.e. 
i?Th=/(cL, ct)/w d . The values of the two quantities are 
very similar in all the cases investigated, the maximum 
deviation, observed once more for the cases of silica 
and selenium, being within a factor two. The column 
Err = \R Exp - RTh\/({REx P + Rrh)/2) reports the 
percentage of deviation between the measured and pre- 
dicted values. Overall, the agreement is very good, also 
in view of the approximation made in deriving Eq. ([5|) 
and, more important, of the experimental uncertainties 
associated with the determination of A and B. 

In conclusion, we applied a recent theory of vibra- 
tional excitations in a disordered elastic medium [29] to 
the determination of the dynamic structure factor. This 
theory, that successfully explained the plateau in the T- 
dependence of the thermal conductivity in glasses [29[ , is 
now found to explain one of the most intriguing features 
in the collective high-frequency vibrational dynamics of 
disordered systems: the k 2 dependence of the Brillouin 
linewidth, i. e. of the acoustic attenuation. The the- 
ory also predicts the existence of a relation between the 
excess DOS (the boson peak intensity) and the sound at- 
tenuation coefficient (the Brillouin peak width). A test 
of this prediction using existing experimental and simula- 
tion data gives an excellent outcome, thus giving further 



4 



support to the validity of the theory itself. Finally, the dots in the insert of [5]) takes place at a frequency much 
theory also indicates that the transition between the fc 4 lower than the boson peak frequency which, in turn, is 
law at very low frequency and the k 2 law for Tk @ (open much smaller than the Ioffe-Regel limit frequency. 
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